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Abstract 

We have studied a model which has been proposed as a regularisation for four 
dimensional quantum gravity. The partition function is constructed by performing 
a weighted sum over ah triangulations of the four sphere. Using numerical simula- 
tion we find that the number of such triangulations containing V simplices grows 
faster than exponentially with V. This property ensures that the model has no 
thermodynamic limit. 
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Introduction 



In the last few years there has been considerable interest generated in a model for quantum 
gravity in which the functional integral over metrics (ill-defined in the continuum) is 
replaced by a discrete sum over random triangulations. The initial proposal |l|, ^ ^ arose 
as a natural generalisation of random surface theories in two dimensions. The results of 
these numerical studies were encouraging and were confirmed by other groups The 
most exciting possibility was the observation of a possible phase transition for a critical 
value of the bare Newton constant. The hope was that a nonperturbative quantum theory 
for gravity could be recovered in the vicinity of this new fixed point. These observations 
were rendered more quantitative by the recent work in which a serious finite size scaling 
study was performed. 

The model is defined from the partition function. 

T{S4) 

The sum is restricted to run over all simplicial manifolds (triangulations) with the topol- 
ogy of . The first term in the action A^4 is just the number of four simplices in the 
triangulation T and this allows us to identify the corresponding coupling clS 3i bare 
cosmological constant. The second term depends only on the number of vertices in the 
triangulation A^q and plays the role of the integrated Ricci scalar - the coupling is then 
essentially the inverse bare Newton constant. 

This correspondence is clear classically from the usual Regge expression for the curva- 
ture associated to any triangle rj^-fc with the extra constraint that the four simplices are 
all considered equilateral 

Tijk = 1tx - cos"^ {^"^'^^^ 

Notice that if the volume is bounded, the number of four simplices shared by a given 
triangle nl^ is necessarily also bounded. This automatically ensures that the model is well 
defined at finite volume - it is a dynamical question as to whether the problems associated 
to the unboundedness of the continuum action return on taking the large volume limit. 

As we have remarked, the analogous model in two dimensions has been studied exten- 
sively, see, for example, the review 0]. It seems clear that at least for central charges less 
than unity, the sum over triangulated graphs correctly mimics the continuum functional 
integrals over the metric including the conformal anomaly. In four dimensions it is not 
at all clear that a simple generalisation, such as the one described above, is sufficient to 
explore the space of metrics. However, it constitutes a simple ansatz which may be studied 
using numerical simulation. 

We may rewrite eqn.|l] in the form 

Z = Y.e-^'''^^{N^,K^) (3) 

The partial sum f2(A^4,Ko) counts the number of triangulations (weighted by the Ricci 
term) with volume A^4. The results we discuss here are concerned with the volume depen- 
dence of this entropy function O (A^4, kq). It is helpful at this point to recall the behaviour 
of the equivalent two dimensional model. 
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Two dimensional gravity regulated using dynamical triangulations has a partition anal- 
ogous to eqn. 0. The number of simplices is now just N2 with corresponding cosmological 
constant K2. The coupling plays no role in two dimensions as the number of vertices A^o 
is strictly proportional to the number of two simplices Nq = ^N2 + x if the Euler character 
X is kept constant (for example S*^). If we sum over all two dimensional triangulations with 
fixed volume we arrive at a quantity u {N2) analogous to Q (A'4, 0) for the four dimensional 
theory. 

There are rigorous proofs @, |^ that this quantity is exponentially bounded. 

00 {N2) ~ e'^^^^ (4) 

This property is crucial for the very existence of the partition function. It implies that for 
a sufficiently large bare cosmological constant K2 > 1^2 partition function will be finite. 
The thermodynamic limit A''2 — > cxd is then obtained by tuning K2 towards this critical 
value The mean volume {N2) then behaves as {N2) ~ {k2 — 1^2)'^ ■ number of 

triangulations were to increase faster than exponentially, it would be impossible to tune 
the bare cosmological constant to approach the large volume limit in a regular fashion - 
the partition function would be dominated by infinite volume triangulations independent 
of the bare lattice parameters. Constructing a continuum limit would then be impossible. 

Thus, it is absolutely essential for the very existence of these higher dimensional models 
that there be such a bound. Unfortunately, there are no analytic proofs available for 
dimension greater than two. If the topology is not fixed it can be shown that the number of 



triangulations increases factorially with volume even in two dimensions [|T0| . The situation 
is made worse by the lack of any topological classification of three and four dimensional 
manifolds. 

Faced with this we have used numerical simulation to estimate the volume of the 
triangulation space. Whilst the previous studies 0, ||, |[ have claimed evidence for 
an exponential bound we believe the issue is of such paramount importance that a very 
detailed study is required. Indeed, the results we shall present favour a very different 
scenario. 



Method 

For an entropy function that behaves exponentially with volume we have argued that it 
is possible to choose the coupling ^4 to fix the mean volume (A^4). In practice this is a 
difficult fine tuning problem. Even under the assumption of an exponential bound, the 
entropy Q {N^, hq) is of the form 

fi(A^4,«;o) ~ A^4^"''^e"4(''o)^4 (5) 

We have included the leading power law correction parametrised by a (ko). In practice 
the power a is negative, so that the partition function is dominated by small or large 
volumes depending on the sign of A/t4 = K4 — kI (kq). 

This problem has been tackled in a variety of ways. We have followed the approach of 
Migdal et al. and added to the action a small correction term of the form 

= 7 {Ni - Vf (6) 
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Replacing the sums by integrals and forgetting for the moment any power law corrections 
it is now simple to obtain a relation between the mean volume (A^4) and the parameters 
in the action. 

{N,) = ^ (An, + 2jV) (7) 

Thus tuning to yield an average volume V yields a measurement of the coupling kKkq). 
The auxiliary coupling 7 merely controls the magnitude of volume fluctuations. We have 
set 7 = 0.005. The presence of power law (and other subleading) corrections gives k1 (kq) 
a volume dependence kI (kq) = {V, kq). The relation eqn. |^ may be rewritten 

fil{N,,Ko) = K, + 2^{{N,)-V) (8) 

In practice we iterate the above relation during the thermalisation stage of our sim- 
ulation and apply it once more at the end of our run to compute our final estimate for 

In this picture the presence of an exponential bound would be signalled by this critical 
cosmological constant nl (V, kq) having a finite limit for large volumes V. In contrast 
^4 {V, Kq) would increase logarithmically in a model for which Q {N4, kq) grew factorially 
with volume (this just follows from the asymptotic result {x\Y ~ g^^'"^). 

Notice that it is sufficient to prove an exponential bound for a single value of Kq - the 
following inequality guarantees that there will then be a bound for any other kq > 0. 

n{Ni,0)<n {Ni, Ko) < exp (aKoNi) n (N^, 0) (9) 

We have used a Monte Carlo algorithm to sample the triangulation space of the model 



- the details are given in . Our code is written in such a way as to make the dependence 
on dimension d trivial - it enters only as an input parameter to the program. 

We have simulated systems from size V = 500 to V = 32000. Typical runs utilised 
on the order of 4 x 10^ MC sweeps with one sweep corresponding to V trial updates. In 
addition we performed a series of runs for both the two dimensional and three dimensional 
models. The results of these simulations could then be contrasted with the equivalent four 
dimensional data and served as an important test of our code. 



Results 

Fig. U is a plot of the critical cosmological coupling (V, kq) against the logarithm of the 
volume for the two, three and four dimensional models at kq = (To improve clarity we 
plot K2 — 0.5 and Kg — 1.0). Clearly, the presence of an exponential bound emerges very 
clearly in the two dimensional {V, 0) is statistically consistent with a constant 

{00) = 1.1249(6) for volumes V > 2000. 

For three dimensions the situation is rather different. The finite volume dependence 
of K3 {V, 0) is large over the full range of volumes analysed. However, as the plot reveals 
there is no strong evidence of a logarithmic component - indeed the best fit we could 
make to the data corresponds to a convergent power law (the solid line in the figure) 
nl = a + bV^. The fit yields a = 2.01(1), b = -3.2(1) and c = -0.28(1) with a per 
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degree of freedom 2.0. Thus, our data in three dimensions favours a bound. Indeed these 
numbers are consistent with the ones quoted in a previous study by Ambj0rn and Varsted 



|12| who give a = 2.06, b = —3.9 and c = —0.32. Their fit derives from lattice sizes of 
V = 14000 and smaller with lower statistics but it is reassuring to see that we are in pretty 
good agreement. We are currently extending our dimension three runs to larger lattices 
to strengthen our confidence in the three dimensional bound. 

The situation in four dimensions is radically different. Clearly, the data support the 
hypothesis that there is a logarithmic component to the critical volume coupling {V, 0). 
A fit of all the d = 4 data to a simple logarithm = a + blnV results in a value for 
b = 0.0315(3) with a per degree of freedom 2.7 (solid line shown). Converging power 
fits simply fail to describe the data. 

To test this hypothesis further we looked at the situation for non zero kq. Fig. |^ shows 
a plot of ^liy, kq) for kq = 0.0, Kq = 0.5 and Kq = 1.0. The inequality eqn.| implies 
that the coefficient of this logarithm should be universal (independent of kq). The leading 
effect of a non zero value for Kq is simply a renormalisation of any exponential terms in 
r2(A^4,0). This is confirmed by the data in fig. 0. Although the curves start out with 
different gradients their large volume behaviour appears to be the same. 

However the plot makes it also clear that the onset of this asymptotic regime is depen- 
dent on Kq - as Kfj increases the curves start off increasingly flat and the logarithm only 
manifests itself for large volumes. 

We found that very long runs were required to thermalise the four dimensional lattices. 
The initial configurations were created by employing only the node insertion move which 
effectively generates lattices corresponding to large values of kq. For the largest volumes 
we employed, V = 32000, we found that subsequent relaxation times were of the order of 
10^ sweeps. This difficulty of reaching true equilibrium was the main factor in determining 
the largest volumes we could reach. It is perhaps a practical demonstration of the results 
reported in |]l3l in which the algorithmic unrecognisability of four manifolds is shown to 



lead to a lack of a reasonable bound on the number of local moves needed to pass from 
one configuration to another. 

Thus our four dimensional data would indicate that the entropy fl {V, 0) has a leading 
behaviour 

n{v,o) r-. {v\f (10) 

If we fit the kq = 0.0, 0.5, 1.0 data for the three largest volumes by straight lines we 
find consistent estimates for the exponent S. These are 6 = 0.027(1), 0.026(1), 0.025(2) 
respectively. We would then assign our best estimate for 6 as 6 = 0.026(5). 



Outlook 

In summary, we have presented results which are consistent with a leading factorial be- 
haviour for the entropy of triangulations of the four sphere Q{V,0). Specifically, the 
number of triangulations of grows like 

n{V,0) ^ exp{aV){V\Y (11) 
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Furthermore, we estimate the exponent 6 = 0.026(5). This rapid growth renders it impos- 
sible to take the thermodynamic (large volume) limit - the partition function for any Kq 
is dominated by large volumes. This in turn implies there is no continuum limit for the 
model. 

We have argued that the presence of large finite volume effects can obscure this be- 
haviour for large values of the inverse bare Newton constant (kq) on lattices that are 
computationally accessible. It is tempting to speculate that the rather rapid shift of the 
pseudo critical node coupling reported in [Q is further evidence for the lack of a well- 
defined continuum limit. The data presented in p is not inconsistent with a scenario 
in which this pseudo-critical coupling diverges as the mean volume approaches infinity, 
leaving the system in an extremely crumpled, degenerate phase. 

It is important to notice also that the term added to help fine tune the cosmological 
constant is now playing a crucial role in defining the partition function. There is now 
no reason to believe that different methods of doing this are equivalent. 

Clearly an extension of this work (with perhaps a more refined method for computing 
K4 {V, Kq)) to larger volumes and node couplings would help to confirm these conclusions. 
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